Measuring Coherent Motions in the Universe 
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We present new measurements of the coherent motion of galaxies based on observations of the 
large-scale redshift-space distortions seen in the two-dimensional two-point correlation function of 
Luminous Red Galaxies in Data Release Seven of the Sloan Digital Sky Survey. We have developed 
a new methodology for estimating these coherent motions, which is less dependent on the details 
of galaxy bias and of the cosmological model to explain the late-time acceleration of the expansion 
of the Universe. We measure a one-dimensional velocity dispersion of galaxies on large-scales of 
o-„ = 3.0ll|^:4^ h-^ Mpc and (t„ = 3.69t|];4^ /i"^ Mpc at a mean redshift oi z = 0.25 and 0.38 
respectively. These values are fully consistent with predictions for a WMAP7-normalised ACDM 
Universe and inconsistent at confidence of 3.8(7 with a Dvali-Gabadadze-Porrati (DGP) model for 
the Universe. We can convert the units of these measurements to 270^41 km/s and 320^41 
km/s respectively (assuming a ACDM universe), which are lower that expected based on recent low 
redshift {z < 0.2) measurements of the peculiar velocity field (or "bulk flows"). It is difhcult to 
directly compare these measurements as they cover different redshift ranges and different areas of 
the sky. However, one possible cosmological explanation for this discrepancy is that our Galax;y is 
located in unusually over, or under, dense region of the Universe. 

PACS numbers: draft 



I. INTRODUCTION 

The last decade has seen a revolution in cosmology 
with the emergence of a standard model for the Universe 
dominated by a strange substance known as "dark en- 
ergy" , with an effective negative pressure. Since the first 
evidence of dark energy in 1998 [l|, [H, there has been 
substantial observational and theoretical research aimed 
at understanding the true nature of this phenomenon. 
In recent years, many authors have started exploring the 
possibility that dark energy, and the observed accelera- 
tion of the expansion of the Universe, could be the conse- 
quence of an incomplete theory of gravity on cosmological 
scales and may require modifications to Einstein's theory 
of General Relativity. 

One of the most direct methods of testing our assumed 
theory of gravity is to perform consistency checks be- 
tween the geometrical expansion history of the Universe, 
as measured by cosmological probes like Type la Super- 
novae and the Baryon Acoustic Oscillations, and the evo- 
lution of density inhomogeneities in the Universe 0-0] ■ 
The growth of structures in the Universe can be measured 
using a variety of techniques, but a popular and direct 
method involves measuring the coherent peculiar veloci- 
ties of galaxies on large scales (i.e., the motion of galax- 
ies after the cosmological expansion has been removed) 
caused by their infall into large scale overdensities like 
clusters and superclusters of galaxies [1, Q. Tradition- 
ally, the local peculiar velocity field, or "bulk fiow" of 
galaxies, has been estimated using samples of galaxies 
where the peculiar velocity of each galaxy has been de- 
termined using secondary distance indicators (e.g. the 



Tully-Fisher relationship [10|-[l^ or Fundamental Plane 
[13l4la |). Alternatively, the coherent motion, or fiow, of 
galaxies can be statistically estimated from their effect on 
the clustering measurements of large redshift surveys, or 
through the measurement of redshift-space distortions. 

Historically, measurements of the local galaxy pecu- 
liar velocity field have been plagued by systematic un- 
certainties [lal, and small sample sizes, which have made 
their cosmological constraints uncompetitive compared 
to other probes of the Universe. However, a recent 
re-analysis of these earlier peculiar velocity galaxy sur- 
veys [17h19| has now provided a consistent observational 
picture from these data, and finds significant evidence 
for a larger than expected bulk flow in the local Universe 
(within 100 /i~^Mpc). This new analysis is also consis- 
tent with new measurements of the local bulk flow using 
the kinetic Sunyaev-Zeldovich effect [13] of massive clus- 
ters of galaxies out to at least 800 Mpc [U , which 
leads to the intriguing situation that these local measure- 
ments appear to be significantly greater in amplitude, 
and on larger scales, than expected in a concordance, 
WMAP7-normalised ACDM cosmological model. Such 
discrepancies with ACDM may give support to theories 
of modified gravity [23| . 

Given the importance of these local large-scale bulk 
flow measurements, we investigated in a previous paper 
[23 | the likelihood of such large coherent flows using an al- 
ternative methodology based upon measurements of the 
redshift-space distortions seen in the clustering of a sam- 
ple of galaxy clusters selected from the Sloan Digital Sky 
Survey (SDSS) [11]. We detected no statistical evidence 
for the coherent flows of the same size, and scale, as dis- 
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cussed above, but due to the sample size the statistical 
errors were large (see [l^l for details). Therefore, in this 
paper, we re- visit our earlier work using a larger cosmo- 
logical data-set and recent theoretical improvements in 
the modelling of redshift-space distortions. The goal of 
this paper is to statistically study the likelihood of the 
local bulk flow measurements within a large cosmolog- 
ical volume of the Universe, well beyond these nearby 
measurements. 

In section II of this paper, we outline in greater detail 
the fitting techniques used in [2^ to formulate the line- 
of-sight smearing effect due to the coherent bulk flow of 
galaxies. In section III, we explain our data analysis of a 
sample of SDSS Luminous Red Galaxies (LRGs; [231) and 
our measurement of their two-dimensional redshift-space 
correlation function. In section IV, we present our results 
from comparing our LRG measurements to our model 
for redshift-space distortions, and compare our results to 
cosmological predictions and previous observations. We 
conclude in section V. 



II. MODELLING REDSHIFT-SPACE 
DISTORTIONS 

The aim of this work is to statistically measure the 
coherent motion of galaxies on large scales. That is, the 
linear velocity of galaxies excluding both the cosmolog- 
ical Hubble expansion and any smaller scale non-linear 
components. 

The coherent motions of galaxies introduce redshift- 
space distortions - an anisotropic feature - into the mea- 
sured clustering statistics. As originally proposed in 
[3^ . a distant observer should expect a multiplicative 
enhancement of the overdensity field along the line-of- 
sight, compared to the transverse direction, due to such 
coherent peculiar motion of galaxies. This "Kaiser ef- 
fect" , as it is now know, can be seen as a "squashing" 
or flattening of the two-dimensional two-point correla- 
tion function (Cs(f, tt)), where the correlation function is 
decomposed into two vectors; one parallel to the line-of- 
sight (tt) and the other perpendicular to the line~of~sight 
(a). Information about the coherent velocities of galaxies 
can then be extracted from the two-dimensional corre- 
lation function via careful theoretical modelling of these 
redshift-space distortion effects. 

In this paper, we utilise several theoretical improve- 
ments to the original redshift-space distortion work given 
in [soj . Throughout this paper, we refer to the origi- 
nal work of [3^ as the "Kaiser limit" , which formulated 
the coherent motions of galaxies as an additional Doppcr 
shift to the cosmological redshift. The "Kaiser limit" is 
suitable to test theoretical models in the linear region 
of galaxy perturbations (on large scales), but over the 
last decade there have been significant improvements to 
the modelling of these effects. Below, we repeat some 
of the formalism presented in our first paper |2J] to aid 
the reader and provide continuity with additional work 



presented in following subsections. 



A. Representation of the correlation function in 
the Kaiser limit 

The observed power spectrum of density fluctuations 
in redshift-space, Pohik,fi)^ can be written as, 

Pob(fc, /i) = [Ps.s, (k) + 2^L^Ps^e{k) + ^i^Peeik)] 

X G{k,n;a^) , (1) 

where 6g and denote galaxy density and velocity fields 
respectively (where Q = 9/aH and 6 is the divergence of 
the velocity fields). In Eq. [T] the function G denotes 
the additional suppression effect due to coherent mo- 
tions which affects the galaxy-galaxy, velocity-velocity 
and galaxy-velocity power spectra equally, and is a func- 
tion of separation (fc), angle between the galaxies (de- 
noted by ^) and has a given smoothing scale, or velocity 
dispersion (cti,), which is the parameter at the heart of 
our analysis and this paper. 

As discussed in our previous paper, we can decom- 
pose these power spectra into a scale-dependent [D) and 
scale-independent part [g) given by, 

P$$(A:,a) = D$(fc)5|(a), 
Ps,s,{k,a) = Dm{k)gUa), 
Pee(fc,a) = A„(%e(a), (2) 

where the subscript $ denotes the curvature perturbation 
in the Newtonian gauge, 

ds^ ==-(l + 2*)dt2+a2(l + 2$)da;^ (3) 

The growth function gf, is defined as gt = b gs^ where 6 
is the standard linear bias parameter between the den- 
sity of galaxies and the underlying dark matter, 5m- As 
we follow the positive sign conversion, is the growth 
function of —0. 

The shape factor of the perturbed metric power spectra 
D^{k) is defined as 

D^{k)^—-Al{k)Ti{k), (4) 

which is a dimensionless metric power spectra at Ueq (the 
matter-radiation equilibrium epoch), and A^^(fc) is the 
initial fluctuations in the comoving gauge and T,^{k) is 
the transfer function normalised at r$(fc — >■ 0) = 1. The 
primordial shape A^^^ (fc) depends on 715 (the slope of the 
primordial power spectrum), as A^^(fc) — A'^{k/kp)"'^~^, 
where is the amplitude of the initial comoving fluctua- 
tions at the pivot scale, kp = 0.002 Mpc~^. The interme- 
diate shape factor T^{k) depends on aj,„ {ujm = ^mh"^)- 
The shape factor for the matter fluctuations, D^ik), 
which is important for both the galaxy-galaxy and 
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velocity-velocity power spectra in Eq. [2] above, is given 
by the conversion from D^(k) of, 



(5) 



where, assuming c = 1, wc can write Hq 
1/2997 /iMpc"^ 

In the Kaiser limit, we would thus assume P{k) = 
P'™(/c), where P'™(fc) is the linear power spectrum, and 
no smearing of the large scale power spectra due to the 
random motions of galaxies within individual dark mat- 
ter halos, or the "Finger-of-God" (FoG) effect seen on 
small-scales in redshift surveys. Therefore, the observed 
compression of (cr, tt) along the line-of-sight due to co- 
herent infall of galaxies around large-scale structures can 
be written in configuration space as. 



4 4 p 

■^9bge + y5e ) 6('')^2(/i) 
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glUr)V4i^^), 



where Vi (p) is the Legendre polynomial and the spherical 
harmonic moment (r) is given by. 



k^dk 



D^{k)ji{kr) 



(7) 



where ji is a spherical Bessel function. 



B. CMB priors 

The coherent evolution parts in Eq.[2]are not generally 
parametcrised by known standard cosmological parame- 
ters. We thus normalise these growth factors at aeq such 
that. 



g'S>{aeq) = 1, 

d9S,„{aeq) 



geia-eq) 



o?lna 



(8) 



Then, we treat gb and a-s free parameters to be mea- 
sured by redshift-space distortions with CMB priors used 
to determine the shape factor _D„ (fc) at the last scatter- 
ing surface. 

The shape factor Dm{k) depends mainly on the fol- 
lowing set of cosmological parameters, ujm, ns-, As and 
rim- While the first two parameters {iOm and ns) are 
well constrained by the CMB at the last scattering sur- 
face (regardless of theoretical model used to describe the 
late-time acceleration of the expansion of the Universe), 
the last two parameters (Ag and 17^) can not be deter- 
mined by the CMB alone. The primordial amplitude As 
has been measured to an accuracy of less than 5% from 



WMAP7 alone, however this constraint becomes weaker 
due to the unknown reionization history of the Universe 
(tighter constraints are achieved by assuming a stepwise 
parameterisation of the reionization history). For the 
matter content of the Universe, the CMB measures 
not ilm, which depends on the Hubble Constant. 

Therefore, given our CMB priors, Dm{k) can be ex- 
pressed as. 



D^{k) 



O A* 

^^ni -^s 



D*^{k : A*g, n*„) . 



(9) 



where Ag and are specific reference values. We use 
the best fit cosmological parameters of the WMAP7- 
normalised ACDM model in this paper, A*g = 
V2.41 X 10^9 and n*„-^ = 0.264. Then the power spec- 
tra can be re-written as. 



Ps,s,ik,a) = [ gt 
Pee(fc,a) = ( ge 



O 4* 
O A* 

^ -^s 



D*^{k : A*g, 0.*^), 



D*^{k :A*s, n*J. (10) 



(6) Here, wc define new growth functions as. 



A O* 
.96(a) = gb[a)- — 

i'm -^s 
A O* 
5e(a) = 9e[a)-^^ ■ 

i Lrn. c 



(11) 



Then the power spectra become 

Ps^sAk.a)=gt'D:^ik:A*sM*J. 

Peeik, a) = 5e'A*„(fc : ^S, ^^rn), (12) 

in which the power spectra are split into two parts, i.e., 
and are determined through our redshift-space distor- 
tion measurements, and is given by the CMB priors 
and the reference values of A*g and ^l'^. 

Thus, we re- write the correlation function in Eq. [Slin 
terms of our new fitting parameters g^ and 5q, 

6(f7,7r)(a) ^ (^gt^ + ^gtgh + ^9e^^ Co('')^o(m) 



+ ^ge'G{r)v,i^,), 



(13) 



where the spherical harmonic moment (r) is given by. 
k^dk 



er(o 



27r2 



■D:^ik:A*s.n*JUkr) 



(14) 



We measure g^ with a normalised D*^{k : not 
ge paired with Dm{k). However, as Pee = g^Dmik) = 
g^Dm{k ■ Ag,Q*^), Pee does not depend on any nor- 
malisation, i.e.. Pee is invariant under the transforma- 
tion between the * (asterisk) variables and those with- 
out above. Thus, we present measured coherent motions 
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FIG. 1. (Left panel) The theoretical 2D-correlation function 5s(fi •"") in the Kaiser limit (unfilled black contours) and using our modified 
Kaiser effect (filled blue contours). Levels of the contours are ^s{a,n) = (1,0.5,0.2,0.1) from the inner to outer contour. Contours in 
the upper panel arc derived from ACDM model with density bias b = 1.82 at z = 0.25, and contours in the lower panel are derived from 
ACDM model with density bias b = 1.92 at 2; = 0.38. (Right panel) The observed 2D— correlation function ^3(0", tt) from SDSS DR7 LRG 
catalogue with the same contour levels. 



not in terms of Qq, but in terms of velocity dispersion 
derived from Pqq. Hereafter, we drop symbol, i.e. 
'wnasierifced' quantities imply 'asteriked' quantities from 
now. 



C. Modifying the Kaiser effect 

It has been pointed out in (sij that the description of 
^s(cr, tt) on large scales (in linear theory) will be modified 
due to dispersion effects in the tt direction. This effect 
was recently investigated using A^-body simulations and 
was found to be important to include [33 |. Therefore, 
in this paper, we have modified the Kaiser limit using a 
Gaussian velocity dispersion term for G in Eq. 1, which 
we will refer to as the "modified Kaiser effect" through- 
out. This modified model has been shown to work rea- 
sonable well on mock catalogues of a ACDM Universe 
(see Hi-Ill). 

In detail, we can write G as 



(15) 



where (t„ is the 1-D velocity dispersion on large scales, 
and is the parameter of interest in this paper. We can 
write a-u as. 



f 

6^ 



-^ee z)dk, 



(16) 



where we see that has units of h Mpc 
k have units of ( /i^^ Mpc)"^ and /iMpc^^ respectively. 



as and 



This is somewhat confusing given (7„ is called a "velocity 
dispersion" in previous literature, but we retain this ter- 
minology here to remain consistent with this literature. 
However, we note that we can convert the units of a.^ 
to km/s using aH(^a)<jy, but this requires a precise mea- 
surement of the expansion history of the Universe (at the 
redshifts of interest), in addition to the determined Pqq 
from our method. 

We therefore modify the Kaiser effect using G defined 
in Eq.[Tni but keeping P{k) ~ P''"(fc). Then the integra- 
tion of in Eq. [S] should be re-expressed, as G changes 
with varying altering in Eg. 1161 to. 



C(0 



k^dkdfik 



COS (^krfik)'Pi(^nip.) 



Here fik is the cosine of angle between k and the pairwise 
orientation. 

For illustrative purposes, in the left-hand panel of 
Fig. [TJ we present our modelling of ^^(o'l''') using the 
standard Kaiser limit (unfilled black contours) compared 
to the modified Kaiser effect (filled blue contours) dis- 
cussed above. In these examples, we have used a galaxy 
bias of 6 = 1.86 and 1.88 at z = 0.25 (upper panel) and 
0.38 (lower panel) respectively (These values of b were 
chosen to correspond to our estimated bias values from 
Table I to illustrate our point). The additional suppres- 
sion in the modified Kaiser effect can be clearly seen on 
large-scales (in the linear regime) at approximately the 
10 to 20% level. 
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FIG. 2. The normalized distribution of 105,831 SDSS LRG galax- 
ies over-plotted with the random catalogue as a smooth red line 
(1,664,948 points). The LRGs are split into two redshift bins 
at the median redshift of the whole distribution {z = 0.32) as 
shown by the blue, dashed line. The low redshift bin spans the 
range 0.16 < z < 0.32, while the higher redshift bin covers 
0.32 < 2 < 0.47. 



III. DATA ANALYSIS 

A. SDSS LRG sample 

In this paper, we have analysed LRGs taken from Data 
Release Seven (DR7) of the SDSS [35|. In particular, 
our LRG sample contains approximately 100,000 galax- 
ies in the redshift range 0.16 < z < OAT with magnitudes 
—23.2 < Alg < —21.2 The 3-dimensional size of the sur- 
vey is best described by the effective volume [66l | 



Veffik) 



n{r)P{k) 
1 + n{r)P{k) 



(18) 



where n{r) is the comoving number density at every point 
r and P{k) is the power amplitude at wavelength k. 
The SDSS DR7 LRG sample used in this analysis has 
Veff ~ \{h^^ Gpc)'^. The LRGs are more biased than 
main sample galaxies and thus have a higher clustering 
amplitude which boosts the signal-to-noise on their mea- 
surements. It is this large volume, and the fact that 
LRGs arc highly biased tracers of the underlying mass 
distribution, which makes them an excellent population 
for studying the statistics of large scale structure. 

In detail, we use the publicly available LRG data 
from [36l - l39j , that include random catalogues and galaxy 
weights that are essential for accommodating the angu- 
lar mask in such large scale galaxy surveys. We have 
converted the observed coordinates of LRGs (RA, Dec, 
redshift) into comoving cartesian distance coordinates 
{Dx, Dy, Dz) using the WMAP7 mean cosmological 
parameters based upon a ACDM model. Therefore, 

,i(r) in Eq. [17] can be re-expressed using a mapping 
from the fiducial WMAP7 normahsed ACDM model of 
/fcgd = kD^/Dl^ (where = (D^ £f-i)i/3)^ 

We note that by assuming a cosmological model for 



the conversion of redshifts and angles to distances could 
introduce an additional anisotropic signal in the corre- 
lation function i.e. the "Alcock-Paczynski" effect (40l |. 
The size of this effect is smaller than the Kaiser effect 
on the scales of interest herein, e.g., if we vary our back- 
ground cosmology, within reasonable limits allowed by 
the WMAP analysis, we see no significant difference in 
our correlation functions. However, we note that a more 
self-consistent approach will be required in the future 
(see [ilHii). 

In Fig. [21 we show the co-moving number density of 
LRGs in our whole sample (compared to the random 
catalogue) and we split this sample into two redshift 
bins as indicated (Bin 1: 0.16 < z < 0.32 and Bin 2: 
0.32 < z < 0.47). We do not split the sample further 
as we wish to keep the thickness of the redshift shells 
{5z ^ 0.15) larger than the scales of interest in our analy- 
sis, thus preventing any aliasing in the redshift direction. 

We estimate the correlation function using the "Landy- 
Szalay" estimator [43 |. 

, DD -2DR + RR 
a<^, ^) = ^ , (19) 

where DD is the number of galaxy-galaxy pairs, DR 
the number of galaxy-random pairs, and RR is the num- 
ber of random-random pairs, all separated by a distance 
(T± Act and 7r± Att. All pairs are weighted using the min- 
imum variance weighting scheme of (45|. Each galaxy is 
assigned a weight according to, 

1 



1 + ni{z)Pw'' 



(20) 



where ni{z) is the comoving space density at redshift z 
and P„, = 4 • 10^/i~'^Mpc'^, as in 0|. The pairs are also 
corrected for fiber collisions, where two or more galax- 
ies separated by less than 55 arcseconds cannot be ob- 
served simultaneously. This is achieved using a secondary 
galaxy weight, proportional to the number of unobserved 
galaxies within the fiber radius. The pair counting was 
achieved using the parallel kd-tree code, NTROPY |47| . 

In the right-hand panel of Fig. [TJ we present our 
measured ^s(ct, tt). The perpendicular component, ct, 
is binned logarithmically into ten bins from l/i^^Mpc 
to 60/i~^Mpc, and the parallel component, tt, is also 
binned logarithmically into ten bins, but from \h^^ Mpc 
to 30ft.-iMpc. 

We use the jack-knife method [l^ to estimate statisti- 
cal errors on ^^(ct, tt), which involved dividing the survey 
into N sub-sections with approximate equal area (and 
thus volume) and then computing the mean and vari- 
ance of ^s(ct, tt) from these N measurements of the cor- 
relation function with the i*'' region removed each time 
(where i = l...A^). In our analysis, we use N = 100 jack- 
knifed samples which were created using an adaption of 
the HEALPix, equal-area, pixelisation codejioj. 

The jack-knife measurements provide an estimate of 
the covariancc matrix for all bins in our correlation func- 
tion. However, this matrix can be noisy, given the num- 
ber of jack-knife samples used, although previous work 
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z = 0.25 




z = 0.38 




20 40 

Mode Number i 

FIG. 3. Eigenvalues (A^) at both redshift z = 0.25 (upper panel) 
and z = 0.38 (lower panel). Dotted lines represent our threshold 
limits. 



[50|,|51| has shown that the number of jack-knife samples 
we use should provide a stable estimate of the covari- 
ance matrix. In order to reduce the statistical noise in 
our covariance matrix further, we follow the same pro- 
cedure as described in [24i and "clean" our covariance 
matrix before inversion j25l427| . In Figure [31 we show 
the distribution of eigenvalues for the decomposition of 
our two covariance matrices and truncate them when the 
eigenvalues decrease below Ai < 2, which avoids any sin- 
gularity in the inversion. This threshold in \i was em- 
pirically determined and, reassuringly, gives reduced 
values close to unity (see below). 



B. Marginalisation over cosmological parameters 

The shape of the power spectra of density fluctuations 
is determined before the epoch of matter-radiation equal- 
ity. Under the paradigm of inflationary theory, initial 
fluctuations are stretched outside the horizon at different 
epochs which generates the tilt in the power spectrum. 
The predicted initial tilting is parameterised as a spectral 
index (ns) which is just the shape dependence due to the 
initial conditions. When the initial fluctuations reach the 
coherent evolution epoch after matter-radiation equality, 
they experience a scale-dependent shift from the moment 
they re-enter the horizon to the equality epoch. Gravi- 
tational instability is governed by the interplay between 
radiative pressure resistance and gravitational infall. The 
different duration of modes during this period result in 
a secondary shape dependence on the power spectrum. 
This shape dependence is determined by the ratio be- 



tween matter and radiation energy densities and sets the 
location of the matter-radiation equality in the time coor- 
dinate. Therefore, the shape-dependent part of is 
well determined, in a model independent way, when both 
ns and LOm are determined at the last scattering surface. 
We thus marginalise our constraints over the growth pa- 
rameters with [ns^LOm) using WMAP7 priors on ns and 
Wm, ns = 0.963 ± 0.014 and (jJ„ 



n 1 QQ/I +0.0056 

U.100^_Q^gQ55. 



C. Non— linear suppression cut— off 

On small scales, there are non-linear effects that are 
not easily modelled via analytical methods, e.g., the in- 
ternal velocities of galaxies within clusters which cause 
the famous "Fingers-of-God" effect seen in galaxy red- 
shift surveys [H^. We investigate here the effect of such 
non-linearities and outline our strategy for mitigating 
them in our analysis. 

First, the contribution of velocity-velocity correlation 
is decomposed from the observed correlation function us- 
ing the anisotropic feature of redshift distorted galaxy 
maps of which information is most extractable at the 
quadrupole moment. Since the leading order of the 
quadrupole moment is the density-velocity correlation, 
what we measure is the cross-correlation between the 
density and velocity fields not the velocity-velocity com- 
ponent. The velocity-velocity correlation is the lead- 
ing order at the hexadecapole moment, which is too 
weak to provide tight constraints on the coherent veloc- 
ity flows. However, both fields (density and velocity) are 
perfectly correlated in the linear regime at least at scales 
k < 0.1/iMpc-^ 

Deviations from a perfect cross-correlation coefficient 
e = 1 (e = PgQ / ^ PggPQQ) uot caused by cos- 
mological effects but rather by non-linear physics. If 
we analyse our data only in the linear regime (where 
e = 1), then gi, and f/e can be simultaneously deter- 
mined from just the density-density and density-velocity 
components. Therefore, we impose a cut-off scale on 
r = Vc^ -|- TT^ at the scale where e docs not significantly 
depart from unity. This corresponds to A: < 0.1 /iMpc~^, 
or Tcut = 15 h^^ Mpc, as seen in N-body simulations [53j . 
We use k ^ 1/r, instead of fc ~ 27r/r, given by the peak 
contribution of convolution in Eq. [T3] due to the Bessel 
function in the integrand. 

Secondly, we also discard more bins along the tt di- 
rection as they are more contaminated by the "Finger- 
of-God" effect. In our model, we use a linear power 
spectrum instead of a non-linear power spectrum as 
in [31I, m, . This should be fine for the accuracy re- 
quired in this paper as the most important correction is 
due to the dispersion effect, which we have included here, 
and also the effect of non-linearities in the power spec- 
tra is only a few percent at fc < 0.2/iMpc~^ (32l. Is^]. 
To estimate the impact of the non-linearity on the two- 
point correction function in redshift space, we use the 
HALOFIT model (Hsj for the non-linear power spectrum 
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FIG. 4. Constraints on gh and av Unfilled black contours denote the result with Kaiser limit, while the filled blue contours 
are the results including the velocity dispersion effect discussed in the text. Contours are 68% and 95% errors. 



with a CDM transfer function [56|. We apply the same 
non-linear power spectrum for the velocity power spec- 
trum, although this obviously overestimates the velocity 
power (see [3J] for a non-linear fitting function of the 
velocity power spectrum). However, we use this non- 
linear model just for discarding small-scale data which 
may be significantly affected by non-linearities, and our 
rejection will be conservative. We find that the effect of 
non-linearity is over 20% at the scale of cr < 5 /i'^ Mpc 
in the redshift range of our data, and discard these mea- 
surements from our analysis. Shown in the right panel of 
Fig.[Tl measured ^s(o', tt) at bins of > 5/i^^Mpc is con- 
sistent with theoretical predictions (compared with the 
left panel of Fig.[T|), but ^s(c, tt) at bins of it < bh^^ Mpc 
does not match well with the theoretical predictions from 
the modified Kaiser effect. In this paper, we thus impose 
a cut-off of (T > 5 ft.""'^ Mpc. 

In summary, we impose two cut-offs in scale to con- 
trol the contamination of small-scale non-linearities on 
our larger scale correlation function measurements. We 
achieve this by not including bins with scales fcut < 
15 Mpc and a <bh~^ Mpc in our fitting procedure. 



IV. MEASURED COHERENT MOTIONS 



A. Consequences of the correction to the Kaiser 
limit 



Using the data and method outlined in Sections II and 
III, we measure gi, and simultaneously from the SDSS 



DR7 LRG correlation functions assuming the appropri- 
ate WMAP7 priors as discussed above. As illustrated in 
Fig. [I] there is excellent visual consistency on large scales 
between our predicted (ct, tt) and the observed func- 
tions. In detail, we obtain a reduced of 0.89 and 0.83 
for the z = 0.25 and 0.38 samples for 12 and 9 degrees 
of freedom. We have 16 and 13 eigenmodes after trunca- 
tion, minus four fitting parameter respectively. In Fig. 21 
we present constraints on gt, and cr„ {a^ is converted from 
measured ge using Eq. I16p . In the Kaiser limit, varia- 
tions of the coherent growth function of galaxy density 
fields amplifies mainly the monopole moment of f s (f, tt) , 
and variations of the coherent motions affects mostly the 
anisotropy of ^^(ct, tt). Measured gi, and pe using the 
Kaiser limit are presented as unfilled black contours of 

Fig.m 



The velocity dispersion in the correlation function 
(ct, tt) of the modified Kaiser effect (Eq. [TT]) induces 
an additional suppression to the (cr, tt) of the conven- 
tional Kaiser limit (Eq. ^ . The variation of coherent mo- 
tion leads to monopole suppression as well as anisotropic 
amplification. As coherent motions increase, the corre- 
sponding gi, becomes smaller due to the increasing sup- 
pression by the velocity dispersion effect. Shown as filled 
blue contours of Fig. 21 the best fit value of gb is shifted 
to lower values and the contours are elongated to smaller 
gb at higher ■ While the measured difference of ct„ be- 
tween the Kaiser limit and modified Kaiser effect is less 
than 4%, the monopole shift of gb due to velocity dis- 
persion with the modified Kaiser effect causes measured 
differences of gb from 10% to 15% presented in Table HI 
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z = 0.25 


2 = 0.38 




(ACDM) 


Kaiser limit Modified Kaiser 


(ACDM) 


Kaiser limit Modified Kaiser 


(T„ in Mpc/h 
aHay in km/s 


(0.411) 
(3.28) 
(294) 


0.380Io.o56 0.377Io.o57 
3.04l°:« 3.0ll0;f6 
27211J 27011? 


(0.422) 
(3.38) 
(293) 


0-478Io.o6o 0.462lo.o59 

q 09+0.49 q f;Q+0.47 
'-'■°^-0.48 ''•"^-0.47 

33llg 320III 


/? 

6acdm 




1 qq+U.U44 -, r,r+U.UbU 
J^-'J'J_0.043 ^^-^0-0.062 
r. 9Q+0.044 r, qn+0.047 
"•^''-0.043 '-'■•J'J-0.048 
1 00+0.066 1 OC.+0.089 
^•^°-0.O65 J^-°"-0.093 




1 q9+U.U4V , 10+U.UV2 
-'-■'5^-0.046 -'-■-'-"-O.OSl 
n qf;+0.048 r, qo+0.056 
U.OU_Q U.O»_Q 

9 nQ+0'074 1 8S+0-11 
z.u»_o 073 J^-°"-0.13 



TABLE I. The measured values of several common cosmological parameters, including gs, gt and av which are discussed 
extensively throughout this paper. The quoted errors on the best fit parameters are one sigma, after marginalizing over all 
other parameters. Parameter values given in parenthesis are theoretical predictions based on a WMAP7-normalized ACDM 
model for the Universe. 



B. Cosmological constraints 

We present herein a measurement of the coherent mo- 
tions using (T„, converted from our measured value of 
(shown in Tabic |I| and shape factor Dm given by 
our WMAP7 priors. At z = 0.25, cr„ is measured to be 
3.0llo;lg Mpc (one sigma errors), and at z = 0.38, we 
find cr-u = 3.691q'1y Mpc. These measurements are 
less dependent on the cosmological model for the late- 
time acceleration of the expansion of the Universe (e.g. 
dark energy) as our CMB priors are fixed at a much ear- 
lier epoch in the Universe, e.g. we use value of A*g and 
constrained at the surface of last scattering. Our 
measurements are also independent of galaxy bias which 
is a major advantage compared to other methods to pa- 
rameterise the motions of galaxies from redshift-space 
distortions. It is true that our measured (t„ correlates 
with Qh, but we do not need to know how to separate b 
and gs^ for determining the coherent motions. 

If we assume H{a) is well-measured, then we can con- 
vert the units of (T„ using aHay (assuming ACDM) to 
give 270115 km/s and 320lll km/s at the two redshifts 
listed in Table ID Although we would not recommend us- 
ing aH(7y for cosmological constraints (without a reliable 
H measurement to an accuracy of a few percent), this 
conversion is more intuitive when discussing coherent mo- 
tions and allows us to compare our values with other pe- 
culiar velocity measurements, e.g., the "bulk flow" mea- 
surements discussed in Section I. 

Traditionally, coherent motions are estimated using 
the /3 parameter d^-dH, 

P,^{k,fi) = {l+/3fi')Ps^sAk)- (21) 

This parameter is equivalent to f/e / 36 in our notation, 
and measured to be /3 = 0.30l!^:5j4^ and 0.39lJ^:^^^ at 
z = 0.25 and 0.38 respectively. As we have measured, 
gt is significantly sensitive to the different redshift-space 
distortion models, and /3 is found to be /? = 0.291qq13 
and 0.36lQ Qlg for the Kaiser limit. Thus, when measur- 
ing /? from the SDSS DR7 LRG samples for use as cos- 
mological constraints, it should be stated clearly which 
model for redshift-space distortions is assumed. 

If we assume a ACDM model, then the galaxy bias can 
be estimated from our measured g^ values. We present 



values for b in Table HI which are 6acdm = l-86lg;Qg3 
and 1.88lo;l3 at z = 0.25 and 0.38 respectively. It is 
interesting to compare our measurements with &acdm = 
1.86 ± 0.07 as measured by [63| using the higher-order 
correlation function of LRGs from the SDSS. This con- 
sistency supports our modified Kaiser formulation, as the 
estimated &acdm using the original Kaiser limit is differ- 
ent by more than one sigma. 



C. Tracing the history of coherent motions of 
galaxies 

In Fig. [51 we present our measurements of ay, the co- 
herent motions of galaxies in redshift space, versus red- 
shift and compared to theoretical predictions. As can 
be seen, our measurements are fully consistent with the 
WMAP7-normalised ACDM model. This represents a 
"clean" test of these cosmological models free from con- 
tamination by non-linear physics on small-scales and un- 
certainties of the galaxy bias determination. These data 
probe the growth history of fluctuations in the Universe 
and are thus complementary to the geometrical probes 
of the Universe. 

In Fig [SI we also present predicted curves for the coher- 
ent motion of galaxies for a variety of dark energy models 
with varying constant equation of state from w = —1.4 to 
—0.6 (dotted curves). All other cosmological parameters 
were fixed to be the same as the WMAP7 best fit ACDM 
values, except w. Using these curves, we can approxi- 
mately estimate that we have a constraint of a{w) 0.2, 
around a mean value of w ~ —1, from our ay measure- 
ments. This is consistent and complementary to simi- 
lar constraints from geometrical observations of the Uni- 
verse [63j . 

Finally, we also provide in Fig [5] a prediction for 
the Dvali-Gabadadze-Porrati (DGP) self-accelerating 
braneworld scenario [gj] . As can be seen, this cosmologi- 
cal model provides a poor description of our observations 
(excluded at confidence of 3.8ct for both redshift bins), 
which support other observational constraints on this 
model from the cosmic microwave background (CMB) 
anisotropy, supernovae and Hubble constant data [65j . 
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FIG. 5. We present the measured evolution of coherent galaxy 
motions at the mean z = 0.25 and 0.38 which are consistent with 
predictions for a WMAP7-normalised ACDM model (solid curve). 
The plotted errors are one sigma as shown in Table I. Dotted curves 
represent predictions from dark energy models with constant w = 
—1.4, —1.2, —0.8 and —0.6 from top to bottom, and dash curve is 
for DGP model. 



V. DISCUSSION AND CONCLUSION 

We present in this paper new measurements of the co- 
herent motions of galaxies on large-scales derived from 
a new methodology for analysing the redshift -space dis- 
tortions seen in the observed two-dimensional two-point 
correlation function of Luminous Red Galaxies from the 
SDSS DR7 sample. Our new methodology is based on 
measuring the scale-independent growth functions, gh 
(galaxy density) and (velocity density), which do not 
depend upon the physics of the late universe (e.g. dark 
energy). They do depend on our knowledge of the early 
Universe, but this can be constrained using information 
from the Cosmic Microwave Background. 

We have determined values of ge from the redshift- 
space distortions seen in the SDSS DR7 LRG data and, 
converting these values into the 1-D velocity dispersion 
(T„, we find = 3.01^0 46 ''^"^ at a mean redshift 
of z = 0.25 and cr,„ = 3.69t^;4^ /i^i Mpc at 2 = 0.38. 
These values for CTu are fully consistent with a WMAP7- 
normalized ACDM model with w ~ — 1±0.2 as illustrated 
in Fig [5l Our observations are however, inconsistent with 
a DGP model for the Universe to high statistical signif- 
icance (> 5cr). Our results provide a competitive, and 
complementary, constraint on these cosmological models 
compared to the usual geometric probes of the Universe. 

We have converted our measured values of into ve- 
locity units (as opposed to lengths presented above) and 



find 270^4° km/s and 320^41 km/s at a mean redshift 
of z = 0.25 and 0.38 respectively, assuming a ACDM 
Universe. As expected, these coherent motions (or ve- 
locity dispersions) are fully consistent with expectations 
from a ACDM Universe. These estimates are however, 
inconsistent with local measurements of the peculiar ve- 
locity field (or "bulk flows") which have recently been 
measured to be greater than these velocities and expec- 
tations from ACDM [13, [3- If the amplitude of these 
local, observed bulk flows were converted (using aHay 
and assuming a flat ACDM model) to the redshift range 
studied here (0.16 < z < 0.47), then we might expect to 
see larger coherent motions. 

It is difficult to perform a direct comparison of these 
different velocity measurements because of the different 
methods and redshift intervals used. For exaniple, the 
lower redshift measurements of [O, [3 or [2l|, only 
probe the velocity field out to 150/i~^Mpc using their 
COMPOSITE sample, while most of the data is within 
a sphere of radius of ~ 60/i~^Mpc depending on the 
weighting scheme used (see Fig 2 of [l3l)- In contrast, 
our statistical estimate of coherent motions on the scales 
up to 60/i~^Mpc are derived from within a large vol- 
ume of the Universe at higher redshift, e.g., 0.5h~^Gpc^ 
(0.16 < z < 0.32) and l.l/i-^Gpc^ (0.32 < z < 0.47) 
respectively, assuming a flat WMAP7 cosmology. There- 
fore, our measurements have averaged over many hun- 
dreds of subregions o f app roximately the same size as 
the volume used by [ITI . llSj (assuming a subregion of 
radius of 60/i~^Mpc). 

Furthermore, there are potential differences in the pos- 
sible directions of the velocity measurements being com- 
pared. The direction of the local velocity measurements 
of [l3, point towards ~ 158,-51 degrees of Right 
Ascension (RA) and Declination respectively on the sky, 
while our SDSS DR7 data is centred at ^ 180, +30 
(Equatorial coordinates) and concentrated in the north- 
ern hemisphere, i.e. our patch is not alon g th e direction 
relevant to the patches used in [13, Ell or |2ll. [22j. A di- 
rect comparison of the directions of the various velocity 
measurements is hard as again our SDSS measurement is 
a statistical average over many subregions of space and 
thus has no directional information. That said, our sta- 
tistical measurement of coherent velocities within such 
sized volumes of the universe is smaller, to high statisti- 
cal significance, than that measured locally around our 
Galaxy. 

A more interesting comparison would be with the mea- 
surements of [m, [121 who find a "bulk flow" of X-ray 
clusters, with respect to the CMB, of ~ 1000 km/s out 
to ~ 800/iYQ^Mpc (or 560/ij~gQMpc for comparison herein) 
in approximately the same direction as discussed above 
for the local measurements [13, Ell- Again, it is hard to 
make a direction comparison in terms of the bulk velocity 
amplitude and direction as our statistical measurements 
are derived from higher redshift [13] , they note that their 
bulk flow velocities peak at z < 0.16 with the possibility 
of higher redshift clusters providing little to their dipole 
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measurements. This is therefore, below our low redshift 
LRG sample data, and averaged over a larger volume, 
than used by [111, [23. 

One possible explanation for the differences we are see- 
ing is that our Galaxy is located in an unusual part of the 
Universe, e.g., in a highly over, or underdense region of 
the Universe. Again, our measurements are obtained at 
z > 0.16, beyond the redshift limits of all these local mea- 
surements [3 HH HI] • Moreover, we can not exclude 
the presence of a large constant large-scale "dark flow" 
across the volume surveyed by our DR7 data i.e., a ve- 
locity dipole with no variation across ~ lh~^Gpc^ of the 
Universe. Such a constant flow would be undetectable by 
our method as the correlation function in redshift space 
is distorted by the divergence of the peculiar velocity 
field, while measuring the bulk flows via the SZ effect 
can include a possible global flow. However, this expla- 
nation would require the "dark flow" of [l^] to extend out 
beyond the SDSS LRG sample (z ~ 0.5) to leave it un- 
detectable, i.e., over a lOOO/i^^ Mpcs in scale. This can 
be tested using the SDSS-III Baryon Oscillation Spectro- 
scopic Survey (BOSS), which will provide redshift-space 
distortion measurements both at higher redshift and over 
larger volumes of the Universe [67| . 

In the future we plan to extend our measurements to 
higher redshifts, in order to extend the history of coher- 
ent motions (e.g. Fig. [S|). As shown in our paper, these 
motions can be measured in an independent way, free of 



some of the problems associated with other measures of 
the growth history of the Universe, and the assumed cos- 
mological model for the late-time Universe, e.g., we are 
able to test not only the conventional dark energy model, 
but also "interacting" and "clustered" dark energy mod- 
els, not to mention the general class of modified gravity 
theories. Unlike other approaches, our measurements are 
free from any possible violation of the consistent equa- 
tions. 
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